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1. INTRODUCTION 


THE present memoir concerns itself with a problem of fundamental import- 
ance in the physics of the solid state, namely the determination of the nature 
of the atomic movements which constitute the thermal agitation in a crystal. 
It would be no exaggeration to say that on our taking the correct view of 
the nature of these movements depends the possibility of our understanding 
the physical behaviour of crystals aright, not only in the thermal field but 
also in diverse other fields, as for example their spectroscopic behaviour 
and the phenomena of X-ray diffraction in crystals. The problem is consi- 
dered in this paper from two distinct points of view. The first method of 
approach is to examine the basic experimental facts in the light of funda- 
mental physical principles and to draw the inferences following logically 
therefrom. The second method of approach is the analytical investigation 
of the problem on the basis of classical mechanics and the general principles 
of thermodynamics and the quantum theory. The result in either case is 
to show that the ideas underlying the well-known specific heat theory of 
Debye and the lattice dynamics of Born and his collaborators are untenable. 
The picture of the thermal agitation in a crystal which emerges from the 
present investigation is fundamentally different from that postulated in the 
theories of Debye and Born. 


2. SOME GENERAL CONSIDERATIONS 


The identification on a quantitative basis of the thermal energy of a 
crystal with the energy of atomic oscillations rests on two fundamental ideas; 
(a) that the oscillators can be enumerated, their total number being three 
times the number of atoms comprised in the crystal; (b) that the energy of 
the oscillators is quantized and for any particular frequency has an average 
value which can be calculated with the aid of Boltzmann’s theorem. We 
proceed to examine these ideas critically with a view to determine their pre- 
cise significance as well as their implications. 


The first of the two ideas referred to above is based on the fundamental 
theorem in classical mechanics which states that all the possible modes of 
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small vibration of a system of connected particles are superpositions of a 
set of normal modes, the number of such modes being the same as the number 
of degrees of dynamic freedom of the system, and further that in each such 
normal mode all the particles of the system have the same frequency of vibra- 
tion and oscillate in the same or opposite phases. In the present context, 
therefore, the term oscillator must be considered as referring to a set of atoms 
vibrating in a normal mode, all of them having the same frequency and the 
same or opposite phases of vibration. 


The second idea, namely that of averaging the energy of an oscillator, 
arises because of the chaotic nature of thermal agitation which involves 
energy fluctuations whose magnitude is determined by the frequency of the 
vibration. The theorem of Boltzmann which enables the average to be 
calculated is a description in statistical form of the behaviour of the individual 
units in an assembly consisting of an immense number of such units which 
are physically identical but differ in their energy states and being in a posi- 
tion to exchange energy with each other form a system in thermodynamic 
equilibrium. In the present context, the units are the oscillators defined 
in the preceding paragraph. Hence, if the application of Boltzmann’s theorem 
is to be legitimate, it is necessary to postulate that in the crystal is present an 
assembly of this nature, viz., a great number of individual oscillators all of 
which vibrate with the same frequency but are in the different energy states 
permitted by the quantum theory. For each separate frequency of oscilla- 
tion, a fresh assembly of that nature is needed to enable the energy for each 
individual frequency to be separately averaged. 


The foregoing theoretical picture of the thermal agitation agrees com- 
pletely with the actual physical picture in the case under consideration. For, 
every crystal is itself an assembly of an immense number of similar and simi- 
larly situated groups of atoms, the internal vibrations of which can occur 
with one or another of the characteristic frequencies of the group. The 
energy of vibration of a group with any one of these frequencies may have 
any of the series of values allowed by the quantum theory and hence would 
vary from group to group in an entirely unpredictable manner within the 
crystal. Such a picture would be in perfect harmony with the basic notions 
of thermodynamics, but would be wholly irreconcilable with any theory 
which seeks to identify the energy of thermal agitation with the energy of 
regular wave propagation in the crystal. But this latter idea is precisely what 
forms the basis of the specific heat theories of Debye and Max Born. It 
follows that these theories are fundamentally misconceived and must there- 
fore fail, together with all their consequences. 
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3. THE NORMAL MODES OF VIBRATION 


The normal modes of vibration of the atomic groupings in the crystal 
referred to above may be determined and enumerated without any difficulty. 
It is obviously not permissible in this connection to regard the unit cell in 
the crystal structure as isolated from the rest of the crystal. For, any oscil- 
lation which occurs in a particular cell will necessarily communicate itself 
immediately to the surrounding cells and set the latter in vibration with the 
same frequency. Hence, to ascertain the normal modes of vibration, we 
have to consider vibrations extending over domains whose dimensions are 
at least as large as the range of the inter-atomic forces. Further, to ensure 
that the oscillations persist unchanged in character, it is necessary to assume 
that the oscillation in each of the unit cells included within such a domain 
is similar and similarly related to the oscillations in the cells surrounding it. 


Remembering that ‘n a normal mode the oscillations of all the particles 
have the same frequency and the same or opposite phase, the foregoing re- 
quirements lead us to the following description of the normal modes of vibra- 
tion: In any normal mode, the oscillations of equivalent atoms in adjoining 
cells have the same amplitude of vibration but their phases are either the 
same or the opposite in the cells adjacent to one another along one, two or 
all three of the axes of the lattice structure. Thus there are 2 x 2 x 2 or 
8 possible situations. As a consequence of the identity of the amplitudes 
of vibration of equivalent atoms in all the cells which need to be considered, 
the 3 p equations of motion of the p atoms in any particular cell contain only 
3p unknowns. Accordingly, it becomes possible to solve the equations 
completely and determine the 3 p normal modes and the corresponding fre- 
quencies of vibration. As there are eight different situations in respect of 
the phases of vibration, the equations of motion would necessarily be differ- 
ent in each of them. Accordingly we have 3 p x 8 or 24 p normal modes 
and frequencies of vibration. 


An alternative and convenient way of describing the normal modes is 
to consider a supercell of the lattice whose dimensions are twice as large in 
each direction as the unit cell. Such a supercell would contain 8 p atoms 
and the total number of degrees of freedom associated therewith would be 
24 p, and this is the same as the number of normal modes permissible. We 
may therefore describe these modes as the vibrations of the atoms comprised 
in the supercell. When this description is adopted, it becomes immediately 
evident that three of these modes would represent simple translations of the 
supercell and would therefore have to be excluded in enumerating the inter- 
nal vibrations of the supercell. We have then left only (24 p — 3) normal 


modes of vibration properly so-called. 
Ala 
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4. WAVE PROPAGATION IN CRYSTAL LATTICES 


As is well known, the classical theory of elasticity leads to the conclu- 
sion that three types of waves can be propagated in any given direction within 
crystal with the velocities determined by the density of the material and certain 
linear combinations of its elastic constants. The analysis which leads to this 
result regards the volume elements of the crystal as simple mass particles whose 
equations of motion can be written down in terms of the tractive forces aris- 
ing from the elastic strains. Hence the theory and its consequences would 
cease to be applicable when these assumptions are invalid and the position 
would then have to be considered afresh from the atomistic standpoint. 
A formal investigation on the latter basis shows that only in the limiting case 
of very low frequencies and of correspondingly long waves are the results of 
the classical theory of elasticity valid. When we enter the region of the 
higher frequencies, the phase velocity and the group velocity diverge from 
each other and the very significant result emerges that the group velocity 
vanishes for (24 p — 3) modes of vibration of the lattice. These modes are 
found to be identical with the (24 p — 3) normal modes referred to in the fore- 
going section. 


To avoid breaking the thread of the argument, the detailed dynamical 
investigation which has been worked out by Dr. K. S. Viswanathan is printed 


as an appendix to this memoir. It will therefore suffice here to comment 
on the significance of the results there derived by him. As is very well under- 
stood, it is the group velocity which is physically significant in all consider- 
ations regarding wave-propagation, since it is a measure of the rate of energy 
transport in the medium. The vanishing of the group velocity for the whole 
series of characteristic frequencies of atomic vibration signifies that there is 
no wave propagation in the real or physical sense of the term when we are 
considering oscillations in which the interatomic displacements play a signi- 
ficant role. It follows that considerations of the kind used in the classical 
theory of elasticity to ascertain the normal modes of vibration of a solid body 
of finite extension are entirely out of place in the atomistic problem. But 
it is precisely such considerations that are employed in the theories of Debye 
and Born to describe and enumerate the movements which they identify 
with the thermal agitation.. Thus, even on the basis of the classical mechanics, 
the ideas underlying those theories are unsustainable and they have therefore 
to be laid aside as being completely unreal. 


5. THE FREQUENCY SPECTRUM OF THE THERMAL AGITATION 


We may sum up what has been said in the foregoing pages by the state- 
ment that the energy of the thermal agitation in a crystal may be identified 





The Nature of the Thermal Agitation in Crystals 167 


with the energy of an immense number of atomic oscillators which together 
fill up the volume of the crystal, the individual oscillators being in the various 
states of excitation corresponding to each of the (24 p—3) frequencies common 
to them all, supplemented by the energy of the low frequency oscillations 
which the three omitted degrees of freedom represent. These (24 p—3) 
frequencies correspond to modes of vibration which are related to the struc- 
ture of the crystal in precisely definable terms. It will be seen that this de- 
scription of the thermal agitation bears no resemblance whatever to that 
envisaged in the theories of Debye and Born; these authors identify the 
energy of the thermal agitation with the energy of waves in immense numbers 
traversing the crystal in all directions and having frequencies which are all 
different from one another, the energy of the wave of any particular frequency 
as well as the sum total of the energy of all the waves taken together being 
distributed uniformly over the volume of the crystal. 


It is noteworthy that in his fundamental paper of 1907 introducing the 
quantum theory of specific heats, Einstein derived his expression for the 
average energy of an oscillator of given frequency assuming, as is indeed 
necessary for applying Boltzmann’s theorem, that the crystal is an assembly 
of an immense number of oscillators all having the same frequency but in 
different energy states and in thermodynamic equilibrium with each other. 
The theories of Debye and Born make use of Einstein’s expression for the 
average energy of an oscillator, while on the other hand the frequencies of 
the waves with which they identify the thermal agitation are all different from 
each other. Thus it will be seen that there is no logical consistency in their 
approach to the specific heat problem. Further, the uniform distribution 
of the thermal energy through the whole volume of the crystal which such 
identification demands is irreconcilable with the fundamental ideas regard- 
ing the nature of the thermal energy of material bodies inherent in the 
Second Law of Thermodynamics, namely, that the distribution of the energy 
over the volume of the crystal should exhibit fluctuations which are the more 
violent the higher the frequency of vibration. 


The foregoing comments do not by any means traverse all the arguments 
and assumptions on which the theories of Debye and Born are based. One 
of the basic objections to their method of approaching the specific heat 
problem is that since wave motions involve progressive changes of phase 
along the direction of propagation and may have any frequency assigned 
to them, they can neither be treated as normal modes nor enumerated. The 
theories of Debye and Born seek to escape this difficulty by postulating that 
the number of wave motions is identical with the number of degrees of free- 
dom of the system, while the choice of wavelengths is determined by still 











168 Sm C. V. RAMAN 


another postulate, e.g., the so-called postulate of the cyclic lattice which is 
claimed to represent the effect of the external boundary of the crystal. Since 
it is obviously impossible to formulate any boundary conditions for the atomic 
movements at the external surface of a crystal, the procedure is clearly arti- 
ficial. But that it is also fallacious becomes evident when we remark that the 
characteristic modes and frequencies of the atoms in the interior of the 
crystal are determined exclusively by the structure of the crystal and by the 
interatomic forces whose range is exceedingly small and hence there exists 
no justification whatever for the assumption which is inherent in the 
postulate that these modes and frequencies are influenced by the presence of 
an external boundary. 


6. THE SPECTROSCOPIC BEHAVIOUR OF CRYSTALS 


It is not proposed in this memoir to discuss how the ideas set forth in 
it impinge on various aspects of crystal physics. We shall however make a 
brief reference to those cases in which facts already well established illus- 
trate the validity of the conclusions reached in it. 


As Einstein emphasized in his fundamental paper of 1907, the theory 
of the specific heats of crystals stands in the closest relation to their spectro- 
scopic behaviour. In seeking to explain the data then available for the case 
of diamond, he made the simplifying assumption that all the oscillators in 
that crystal have a single common frequency. He recognized however, that 
in general a crystal would possess several characteristic frequencies and that 
the expression for its specific heat would involve a summation over them all. 
In the particular case of diamond, the (24 p—3) frequencies indicated by the 
present theory are reduced from 45 to 9 in number by reason of the cubic 
symmetry of the crystal. One of these 9 frequencies represents the triply 
degenerate oscillation of the two lattices of carbon atoms against each other, 
while the other eight represent either the longitudinal or the transverse oscilla- 
tions of the layers of carbon atoms in the cubic and octahedral planes of the 
crystal. In numerous papers published in these Proceedings in recent years, 
it has been shown that all these nine frequencies can be ascertained by 
spectroscopic observation and that the specific heat of diamond can be 
successfully evaluated with the aid of the data thus obtained. 


But the case of diamond does not stand alone in this respect. In recent 
years, the spectroscopic behaviour of numerous crystals have been studied 
by diverse methods, e.g., the scattering of light, absorption in the visible, 
ultra-violet and infra-red regions of the spectrum, and in particular cases, 
also their luminescence. The whole of the evidence which has thus come 
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to hand indicates that the vibration spectrum of a crystal consists of a dis- 
crete set of monochromatic frequencies and that where there is any departure 
from this state of affairs, an adequate explanation is forthcoming. 


7. X-RAy DIFFRACTION IN CRYSTALS 


That the atomic vibrations in crystals appear with precisely defined 
frequencies and with modes bearing a precise geometric relationship to the 
architecture of the crystal has a most important bearing on the subject of 
the diffraction of X-rays by crystals. Here again, a particularly striking 
illustration is furnished by the case of diamond. When a plate of this crystal 
is traversed by a narrow pencil of X-rays and the resulting Laue pattern is 
recorded on a photographic plate, it is found that in addition to the usual 
Laue spots, other sets of sharply defined spots appear on the plate, corres- 
ponding to each of the monochromatic components in the incident X-radia- 
tion. This phenomenon which was discovered in the year 1940 by 
Dr. Nilakantan and the present writer using octahedral cleavage plates of 
diamond has since been thoroughly re-investigated at this Institute by 
Mr. Jayaraman and the writer. The results are found to be in complete 
accord with the hypothesis that the incident X-radiation excites the triply 
degenerate oscillation of the two lattices of carbon atoms in the crystal against 
each other. These movements occur along the directions of the three 
cubic axes and as a consequence, all the lattice planes of the crystal are dis- 
turbed, giving rise to dynamic X-ray reflections by them. Three sharply- 
defined spots corresponding to each monochromatic component in the 
incident X-ray beam are demanded by the theory for the dynamic reflections 
by each of the (111) planes of the crystal, since these planes are equally 
inclined to the cubic axes. These reflections are actually observed in the 
directions indicated by the theory with suitable settings of the crystal. Other 
sets of lattice planes in the crystal also give dynamic reflections in the direc- 
tions demanded by the theory, but since such planes are unequally inclined 
to the three cubic axes which are the directions of movement of the atoms 
in the excited oscillations, they are not all equally comspicuous. 


Here again the case of diamond does not stand by itself, since analogous 
phenomena are also exhibited by other crystals. But it is not necessary 
here to enter into the details of these cases. 


8. SUMMARY 


The principles of thermodynamics and of the quantum theory indicate 
that the thermal energy in a crystal is identifiable with the energy of an 
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immense number of atomic oscillators which have their frequencies of vibra- 
tion in common but are in different states of excitation. Dynamical 
investigation of the problem shows that if the crystal contains p interpene- 
trating Bravais lattices of atoms, there are (24 p—3) frequencies of vibration 
common to all the oscillators. An investigation of the propagation of waves 
in crystal lattices shows that the results of the classical theory of elasticity 
are valid only in the limiting case of very low frequencies and proves further 
that the group velocity of the waves vanishes for (24 p—3) frequencies which 
may accordingly be identified as the characteristic frequencies of atomic 
vibration in the crystal. Since there is no wave-propagation in the real or 
physical sense except in the very lowest range of frequencies, the identifica- 
tion of the thermal energy with the energy of waves traversing the crystal 
in all directions and having frequencies all different from one another which 
forms the basis of the specific heat theory of Debye and of Born’s lattice 
dynamics ceases to be justifiable. Other aspects of these latter theories 
are also examined and it is shown that the assumptions which they involve 
are untenable and that in consequence the theories have no claim to validity. 
The experimental results which confirm that the atomic oscillators have 
specific frequencies in common and modes of vibration specifically related 
to the crystal structure are briefly set out. 
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APPENDIX 


We shall consider the most general case of a crystal whose structure 
includes p interpenetrating Bravais lattices, all of which are geometrically 
similar and similarly situated in the crystal arrangement. An unit cell of 
the crystal will therefore contain p atoms whose masses we shall denote by 
M,, Mg,---+ Mp. 


Since a crystal is triply periodic in space, any unit cell can be reached 
from any other cell of the crystal by the application of three primitive 
translations or integral multiples of them. We denote the three primitive 
translations by the three vectors d,, d, and d,. In the same way, we denote 
the three primitive vectors of the reciprocal lattice by b,, b, and b,. In what 
follows, the symbols s and o will be used to denote particular cells in the 
crystal structure while r and p will indicate particular atoms in the unit cell; 
X, Y, Z are three mutually perpendicular directions chosen as the co- 
ordinate axes along which the displacements of the atoms from their equi- 
librium positions are resolved. We represent the displacement-components 
of any atom, say the rth atom in the sth cell, by means of the symbols 
dzrs> Tyrs> Yzrs- ‘We can now write the equations of motion of the atoms 
in the sth cell and they take the form 

— Mrq xrs = ~ kupe ype C = De tee ) 


r=1,2,....p (1) 


The force constant k*f¥ occurring in (1) expresses the proportionality 
existing between the x-component of the force acting on the atom (r,s) to 


the displacement in the y-direction of the atom (p, ¢) which gives rise to this 
force. 


Since the crystal is composed of p different Bravais lattices of atoms, 
we can reasonably expect a wave of a given wavelength and frequency to 
be propagated with different amplitudes inside the p different lattices. To 
solve the set of equations (1), we therefore assume solutions of the form . 


ders = Agr eot-0-s) ( =X, Y,Z ) 
2ni 


a Dn ex (vt—e.s) 


(2) 
p=], 2, ...2 

which are plane waves of wavelength A and circular frequency w propagating 
in the direction of the vector e of the crystal. If we resolve the wave-vector 
a along the three axes of the reciprocal lattice and write a = 6,b, + 0,b, 
+ 6b, then the equations (2) can alternatively be written as 


Gers = des etlwt—S191—S292—Sa9s) (2’) 
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where 5,, S, and s, are the components of s along the three Bravais axes of 
the crystal. Substituting the equations (2) in (1), we obtain 


x= X, Y,Z 


-_ ars pid.o—-S 
Mw? Agr — ~ Ayp (2 yPo e ) ( = l, y a eee -p (3) 


Equations (3) are a set of 3 p equations in the 3 p amplitudes Az,,'A,,, 
Az, Axs,----Azn of the waves. By eliminating the amplitudes from these 
homogeneous equations, we can obtain a determinental equation containing 
3 p rows and 3 p columns whose elements are the coefficients of the various 
Ayp’s in the equations (3). By expanding this determinant, we can rewrite 
this equation in the form 


Spw®P + syw9P-24.... 2, + Ssp = 0 (4) 


in which Sp, S;, ....Ssp are functions of the coefficients of the amplitudes 
in the equations (3), and hence are also functions of the wave-vector of the 
waves. The roots of the equations (4) are consequently functions of the 
wavelengths of the waves, in other words, the waves are dispersive. Since 
an algebraic equation of degree 3 p has 3p roots, for each value of a we 
obtain 3 p permissible values for w which are the roots of (4). Thus there 
are 3 p waves of a given wavelength moving in any direction of the crystal. 


We shall now consider any one of the 3 p roots of the equation (4). If 


this is denoted by w,;?, then w,? expressed as a function of the wave-vector a 
has the following properties: 


(i) w,? is a periodic function of 6,, 6, and 0, with periods 2 7 for each 
one of these variables; 


(ii) w,” is a real function of 6,, 0, and 4; 

(iii) w;? is an even function of the wave-vector a (i.e.), w_? (a)= wp? (—a). 

Condition (i) follows from the fact that the coefficients of the ampli- 
tudes in (3) are of the form 


> k+l? exp. ia.(o — s) or 
Px kupe exp. i(o, _ $9, ot. copy — $205 ot. 03 — 593 
and hence are periodic]functions of 0,, 6, and 4; with periods 2 for each 


of them. The frequencies of the waves which are functions only of these 
coefficients are therefore also periodic functions of 4,, 8, and 43. 


Condition (ii) can easily be understood from physical considerations. 
Condition (iii) follows from the fact that the frequency of a wave of given 
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wavelength and amplitude should be the same both for waves travelling in 
one and the opposite directions. Alternatively, these two conditions are 
consequences of the hermitian property of the coefficients of the amplitudes 
of the waves in the right-hand side of the equations (3). 


Since w;? is a periodic function of 6,, 9, and 8, we can express it as a 
three-dimensional Fourier-series in the form 


+co 
wr? => zr zx a Ann, et (N191+N292+Ns93) (5) 


mm, 1, n= — 


where the summation extends over all integral values of n,, n, and ng, from 
— co to + ©oo~ 


If we use a single symbol n to denote the triplet (”, n, n,) we can rewrite 
(5) also in the form 


wp? = J Ay e-" (6) 
Since w;,2 is real, the coefficients of e@(™6:+"202+%62) and e-i(nie14n2g2+n2gs) 
in (5) should be complex conjugates. Hence we have 
An = An (7) 
Now by writing — a for a in (6), we get 
w2 (— a) = J An en (8) 


Since from condition (iii), w,? (— a) = w,? (a) we see by a comparison 
of (6) and (8) that 


An = A_n (9) 
Combining (7) and (9) we get 


An = An = An (10) 
or Ay is real. 


We can therefore rewrite (5) as 
we? = Agog + 2X Anynin, COS (19, + MgFg + 1345) 
= Agoo + 2 5 An COS a.n (11) 


where the summation now is over all lattice points on and above any one 
of the lattice planes d, = 0 or d, = 0 or d, = 0. 
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We shall assume that the series obtained by differentiating term by term 
the sum on the right-hand side of (11) converges uniformly for all values of 
6,, 8, and @, in the interval (0, 27). Then, we have 


2wK - =—2zAnensin an (12) 


Now, when 6@,, 9, and 0, each takes one of the values 0 or 7, sin a.n 
(i.e.) sin (”,9, + m9, + n39,) is equal to zero. Hence if 
wK ~O (13) 
for these values of 6,, 8. and 43, we find from (12) that 


-_ =0 whenever 0,=Oor7; 
* . 
6, = Oor7; 
6, = Oor7 (14) 


When 0,=0; 0,=0; 0, =0, it can be proved that three of the roots 
of (4) will become zero and hence the condition (13) is not satisfied. These 
are the limiting cases of elastic vibrations of long wavelength in the crystal 
and for them, since #, > 0 as a>0 


it = da (k = 1, 2, 3), 
and therefore the group and wave velocities become identical with the velo- 
cities of propagation of the elastic vibrations in the crystal lattice. 


The eight set of cases considered above now yield a total number of 
24 p frequencies. Leaving out of these the three degrees of freedom per- 
taining to pure translations and for which, as we have seen, the group velocity 
is non-zero being equal to its wave velocity, we get (24 p—3) frequencies for 
which the group velocity of the waves vanishes. The modes associated with 
these frequencies now possess a simple geometric description. By substi- 
tuting the values of 0 or 7 for 4,, 6, and 43 in (2’), it can easily be seen that in 
all the (24 p—3) cases referred to above, equivalent atoms in adjacent cells 
vibrate with the same amplitude and with the same or opposite phases along 
one, two or all three of the Bravais axes of the crystal. Since the group 
velocity is zero, any disturbance with these frequencies remains centred in 
the region of its origin. Further, the phase relationships for these (24 p—3) 
modes are exactly the same as those demanded for the normal vibrations of 
a connected system of particles. 
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Received August 5, 1955 
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INTRODUCTION 


THE distribution of thermal neutrons in a medium due to time-dependent 
fast neutron sources has been investigated on the age theory approximation 
by Gallone and Salvetti.1 They have also obtained solutions for certain 
specific geometries using the ‘symbolic operator’ formalism. 


In this paper we shall discuss the multigroup approximation for slowing 
down process and obtain the age theory approximation as a particular case 
resulting from the number of groups tending to infinity. It will be shown 
that the resolvant operator used by Gallone and Salvetti is just a particular 
case of a more general resolvant operator, namely the multi group resolvant 
operator. Finally we shall obtain explicit solutions on two group 
and age approximations for an infinite geometry as well as for a finite rect- 
angular parallelopiped. The age theory expressions for a finite geometry 
have already been given by Gallone and Salvetti and are included here for 
completeness. 


THE MULTIGROUP APPROXIMATION 


In the multigroup approximation the neutron energy range from source 
to thermal is divided into a finite number of intervals. The neutrons in 
each interval of energy or ‘group’ obey the ordinary diffusion equation 
and it is assumed that the neutrons in any particular group suffer collision 
without energy degradation till they have existed in that group for a time 
equal to the slowing down time for that energy interval, at the end of which 
they go over to the next group. 


We shall divide the energy interval from source energy (which we shall 
throughout consider monoenergetic) to the thermal energy into (n — 1) 
groups, so that we have nm groups including the thermal group. If Lx be 
the slowing down length, and ¢, the slowing down time of the neutrons in 
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the group k, Ly, the diffusion length and t, the mean capture lifetime of 
thermal neutrons, then we have the following n differential equations for the 
n groups: 


t d 
L;? v7? Pk (r, t) ini Pr, (r, t) + — Pry (r, t) = tk PY; Pk (r, t) (1) 
(k = 1, 2, 3,....,n) 


where px (r, t) is the density of the kth group neutrons at the point r and 
time t.p, (r,¢) is the source term and 4% is taken to be unity. Writing 
Pk (¥, t)/tk = ox (r, t), we get for the thermal neutron distribution the 
equation 


(1-2 2 + te ay) [(1— Lat V? + tn gy) bn] =H 000 
Q) 


From (2) we get as the formal solution of the thermal neutron distribution as 


"T 202 »\" 202 2\7 
én(t, 1) =H (1— Lit V8 + tes) (I Lat V8 +t 3) bol 
(3) 
n—1 \-2 ; 
We define JT (1 — L;? V? + te =) as the time-dependent n-group iesol- 
k=1 
vant operator. 


For the steady state we have instead of (1) the following equations 


Li V? p(t) — pe (®) + gpk (t) = 0 4 


(k = 1, 2, 3,...., 0) 
which yields for the thermal neutron distribution 


TT (1 — Ly? V2) (1 — Ln? V2) $n ®)] = $0) (5) 


ke1 


Hence the formal solution for the thermal neutron distribution can be 
written as 


bn (F) = 1 — Le? V9 — La? V9 bo) 


= [z Ax (1 — Ly V*) |eo (r) (6) 


k=1 
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where A;, is the coefficient of the kth term of the expansion of the continued 
products x (1 — Ly? A*)"* into partial fractions, where A? is an algebraic 
variable. This expansion is possible since the operators (1 — L,? V7?) 
(k= 1, 2,....,) are commutative. Here 7(1 — L;? V*)' is the n-group 
resolvant operator for the steady state. a 

AGE THEORY APPROXIMATION 


In the age theory approximation we assume the energy loss of neutrons 
from source to thermal energy is a continuous process. This is equivalent 
to an infinite group theory. Hence the resolvant operators in the age theory 
approximation is obtained from the n-group resolvant operators by making 
noo. Then we get 


(i) for steady state: 
L (7 a-u2 vy") = en (7) 
aco k=1 
if the sequence XY L;,” tends to the finite limit t+, when k > oo. Here 7, is 
k 


the square of the mean slowing down length from source to thermal energy. 
It is called the ‘age’ of the source neutrons to thermal energy. 


(ii) for time-dependent case: 


n—1 i ‘ Pa) 
ai $29 <e : = e%V*-t.— 
woco Neat [! —— se a jee 8) 


if the sequence Y Li? > 7) as kK 00 and J ty > ty as k —>o0. 
k k 


Here 7» is same as in (7) and ft is the mean slowing down time from source 
to thermal energy. Equation (8) could also be obtained from the age 
equation (see Appendix A). 


The thermal neutron distribution in the age theory then becomes 
$n (r) = eV" (1 — Ln? V7)" oo () (9) 
for the steady state and 


$n (Fr, t) = ereV*-t = (1 —L,? V7 + tn 2" do(r,t) (10) 


for the time-dependent state. 
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If ¢o(r, t) = Do (r) f (0), then 


2 2 : g) 
eV 5 by (Ht) = eV" Dy (n) e a¢ f(t) 
; > t 2 >2 3 3 93 
= eV ®, (r) {1 — ty yt +a sa ay sate} £0 


=e 70,0 [fO-4fO+Pro-frro+-] 


The expression inside the square bracket is the Taylor expansion of 
S(t — t ). 


Therefore, 
sai ae 2 
e7eV ty ye Go(r, th) =e?Y Dy (r) f(t — tf) 


= eV" bo (r, t — tf) (11) 


Hence on the age theory the thermal neutron equivalent source is 
e7*Y" do (r, t — ts) where ty is the slowing down time. In the multigroup 
theory the neutrons in each group diffuse for a finite time equal to the slowing 
down time for that group before a transition to the next group takes place. 
Hence an initial burst of neutrons has a spread in time after each group and 
does not maintain the initial time characteristic. In the age theory these groups 
are infinitely small and hence the neutrons exist in any group only for an 
infinitely short time. So there is no diffusion in any group and hence the 
time characteristic is maintained throughout the fast region. We shall call 
this phenomenon time diffusion. Obviously the best approach is the multi- 
group theory with number of groups somewhere around the number of 
collisions. When the number of collisions is large, this would not differ 
very much from the age theory, but in the case of very light moderators like 
water, the continuous slowing down model is a poor approximation, since 
the number of collisions is small. 


ExpLicir SOLUTIONS FOR THERMAL NEUTRON DISTRIBUTION 


The thermal neutron distribution is given by 


> 
én (t,t) =F (V% 5) do) (12) 
where F ( Vv, =) is the appropriate ‘ diffusion’ operator. The solution 
of (12) is obtained by expanding ¢y (r,t) and ¢(r, 7) into a: series of 
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orthonormal functions A, (r) exp.(ivt) in r and ¢. Here A,(r) are the 
eigenfunctions of the equation 


V7 Ay (r) + w2, A, (F) = 0 (13) 


with the boundary condition that A, (r) and all its spatial derivatives vanish 
on the boundary to which A, (r) is defined. w,* are the corresponding 
eigenvalues. , 


For a discrete eigenspectrum [this is the case when A, (r) is defined within 
a finite boundary] we have 


dn (, t) =I Au (r) [B, (») et dp 
M=0 —co 


co +co 
do (t,t) = 5 Av) f Ay (vy) e”* dv (14) 
H=o0 ——e 
Here A, (v) and B, (v) are the coefficients of the expansion. A, (v) can be 


determined from the nature of the source term ¢)(r, ¢). Substituting (14) 
in (12) we get 


co +00 sia it 
i 2 9 , 
5 A, (r) f B, (») e! dv =F (v3, =) ) A, (r) JA eitdy. 
H=0 -cO poe a= 


(15) 
From the orthonormal properties of A, (r) we have, since 
F(v% 3) uMett =F (— a, i”) AO eh 
B, (») =F (— w,”, iv) A, (») (16) 
Hence 
co +09 
$n (tr, t) = FA, (8) f Ay (») F (— 2,2, iv) e”* dv (17) 
u=0 a 


This completes the formal solution for finite boundaries. 


When the eigenspectrum is continuous, the summation is replaced by 
integration. Let k*? be the continuous eigenspectrum of the equation 


V? A(k, r) + 2A (k, r) = 0. 


Then A(k, r) is a plane wave’ exp. (ik.r). Hence we get 


’ d - ; F sy 
F (V4 5) AG ne” = F (— #, iv) cite ot (18) 
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so that 
B (k, v) = F (— k?, iv) A(k, v) 


Therefore the thermal neutron distribution is given by 
+00 
dn (r, t) = / eik.r dk / A (k, v) F (— k?, iv) e*”* dv (19) 
(k) —co 
This completes the solution for continuous eigenspectrum (i.e., infinite 


boundary). 


We shall now evaluate (17) and (19) on age and twogroup approxima- 
tions for a finite as well as for an infinite geometry. We shall consider an 
instantaneous point source of fast neutrons as well as a point source of 
finite width in time. 


For an instantaneous point source we have 


$o (F, t) = $03 (8) 3 (t) (20) 
where ¢, is the source strength. Then 

A (k, v) =.A (k)/27 = 1/1674 (21) 

A, (vy) =A, /20 (22) 


For a rectangular parallelopiped of sides a, b, c, origin at one of the corners 
and source at (x’, y’, 2’) 


* = _ 8b .. bx’ ... my’ ... nz’ 
A,= Aimn = we 1 | Sin —— sin — 
Nu (r) = Amn (x, ys z) = sin ben sin wn f sin =— 
a b c 
2 m* 2 
w,” = mn" =7 EF + + 7 | (23) 


For a point source of finite width in time we have 


do (tr, t) = $98 (F) (I(t) —1(¢— Ad) (24) 


where I(t) is the Heaveside unit step function. Then 


At 
A, (>) = - § el” dx (25) 


At 
l = 


A(k, ») = jga f ev dx (26) 


Q 
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I. The Age Theory Approximation 
The diffusion operator is given by 


F (v? ss) = envi 2 (1 — Ly? V2+t d\7 
°- 2 " ot 
I.1 Finite Geometry 
(a) Instantaneous source 


co +00 
de 1 oT oy *_ivty evt 
tn (r, t) ae 2 P Ay An (7) f al +; Wy v2 L,? + ivtn) + ist) * 
M=0 -S 


- w “ohn 
=F DUE DA a @ erent oo AR 1) 
In 





&=0 
(27) 
where I(t — ty) is the Heaveside unit step function. 
(b) At source: 
Au (v) e Tow? et ents 
anal’ p2 . of “C+ 0,8L,? + it, 
. f ev (t—ty—x) 
Tow, x x 
a Dy Mee ur, (r) dx fate * 
1 + wy L,? 
2 1 — ge ty) 
= > I(t — ty) A, A, (7) e-7omm Soa 
a ‘hed |) 
(0 <t<AJAd) 


s l+w, 1+wy?L,? 
= ye t4) Ay A, (r) entmn® e- —F—™ (t— (em At -1) 
“ 


(t> At) (28) 
1.2. Infinite Geometry 


(a) Instantaneous source: 





ntno~ Loe te heat 
nam fe fA oat Sim & 


(k) 


: +i sir Mat 
= ch f & core tee wwe 
n 


(k) 
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L,? he t) 
_ etttnite do 1(t — 4) f sie oe [+0 + Pe sin kr 





27? rin Vir 
by e-t-titn T(t — ty) x 
m — tien }- - 
izzy, 4 [ t) + ETH 
n 
2 =e —3/ 2 
x [4 (+ 9) (29) 


(b) At source: 


I ik k* dk Ca @ ely (t*-2) 
= Tr To 
bn ("5 1) = 16 qa f oe iY xf (i + k®L,? + vty)” 





At 
1 ; ; 1+ KL,2 
~ oe f eik.r otk” dk f l(t*— xe 3 “*-») dy 
(k) i) 
1 f ent? otk. [1 (“== 1* ye] 
= 3 27. 2 —, ~ I (t*) dk 
™ aw oo (0<t<AJ)d) 


I 


1 eT oikr 1 + K*L,? 1+ AL,7\ , 
Qr3 if i+k*L,? [e(- te ) Om se 1] e- ( tn ) ' 
ws I (t*) dk (30) 
where t* = t — fy. 


Now 


co 


ek. = g-ak? dk=2 e-ak? k — 
f iro” "Ser s4 i > 


(k) 


-* ke-*** Sin kr An 
fi TSeLS dk = —- I (a, Ln) 


where I (a, Ly) has been shown in the appendix to be 


(ost) aie Lem fier (YE — a) 


vem fier (saz) 


erf (x) = <. f eo du 
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Hence (30) becomes 


bn (tt) = S55) [I(r Ln) — et T(@, Le) | 


0 <t<Ad) 


= EO) [ett 1 (ry Ly) — et TL] >A) 2 


where 
aa <3 sai Wats be hs t) 
nr n 
*=(t—t,); t'* =(t— t7 — Ad). 
II. Two Group Theory 
The diffusion operator is given by 


F (vs, 7) (i-Le We ay ( — Ent? + tn =) 


where L;* = 7», the ‘ age’ and ¢y is the mean slowing down time in the fast 
group. 
II.1. Finite Geometry 

(a) Instantaneous point source: 


400 
a 1 eit 
tn (r,t) = 5. 2 Ay Ay (r) f (+ o2LP+ivt) I-+02L,2+ iit, 4 
mm —oo 


= Ay As (7) Ee eos iT 
id tn tf ®, Tr —@D j 








(33) 


Where 
of = 1 + w,*L;? «+ Bie. 1 + w,?L)? 
a t f . 7 tn 


(b) At source: 


a 1 —_ o et 
$n (1, t)= 20 x wif asfe ia *L in? +ivtn) (I+, *L5?+ivts) 
be 


ev” (t-x) 
m Dy Auda f: of [I-Fo, "Ln? +ivt,) Ue, 2L y+ ivi) 4 
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t 
= AyAy (7) fi (t me x) emu" (t—x) __ e®,! (t-x) 


——_——— dx 
— tn lf - [o," = by!) 


AA, (7) [ l—e Dit 1—e Pt 
- ae ——— — = 1 O<t <a 
2X tn tf @,.” = @,!) Oy! Du” ( AA) 


th tt (Gu” — 5) bi Du" 





_ > AuA, (7) Eee eh: ge Nar a 
& 


(t> Ad) (34) 
where ®,”-J are same in (33) and ¢’ = t — At. 


II.2. Infinite Geometry 


(a) Instantaneous Source: 


+00 . 
_t ( & iia 
n(r,t)= jer ) okt dk f (LAL, +ivty) C-FRALz? + vty © 
(k) = 


1+ WL? , 1 + KL, 
poe! f eik.r [< . Weiss. en =x | dk 
88 ty (1 + k*Ly?) — th (1+k?L;? 


(k) 


1 ae eik.r e-ak? _ eik.r en ank* 
= |< f “Te e™ f Tae dk| 
(k) (k) 


1 
= xaim | et Ley Ve) — Ten Ve) | (35) 
where 
_ Lt aia ig ae 
In p tf 
2 a 2 


(b) At Source: 
ws 


t +co , 
l : etv(t-2) 
A ii tk. k ——+35 9. OOS SS es SC - 
On (r,t) ron J “idtien J & sf (-FML,? +ivty) FRLp pint) @ 


LHL? Gy M+ RL 
oe etkar dk ; - «¢ ( O_o ; a 
= 8nd) tty asi 

( 0 





k) 


(* tet) a ess) dx 
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for 0O<t<AJAt 
[ 1+K'L,? 
e 


1 , . 1 
dn(r, )= ol f tretkr e-tity e-ark (B a WW > - dk 
(i) 


[ 1+K'L,? ora: ] 
ibe tr etker entitn e—ank" e tn ] 
" (B + Ak*) (I + PL a) 4 


(k) 


Pt oar 10, vw“) * ie Sila Pr sar 


— 1, Vun) tr ae , + 10, Vus) ar 

















petits 
— I (an, Vp) tn oH I (a, Vu) Y me 
entitn tits 

+ U(Gns Vin) tm Ge, — Log, ving) HE) (36) 

for t> At 
eiker ty 1+/L,? +L? 
dn (r,t) = ame a > uk?) [7 + Lp ly r( e ty ae 1) 
atts 3 sl 
-TEROIe t(e t, A' —1)| dk 
1 t’itn 

= aang [ te M— an, Va) — roe “Tes vu) 

— ty A on’, Vain) + ty ME ot Cay’, vs 

petit tit, 
— tn Tons Vier) + ty Prat ‘I (ay, Vu) 
+ tn ae ges (an, Ven) — tf a pt e"4T (az, Vf) ] 
, 37 
where ii 
an = Ly" t/ty An’ = Ly? t'/tp 
’=t—At 
= Lf*t/t; as’ = Ly*t'/ty 
Bn, = Ly? B eek 


pe = L;? & = (Ln'ty — Ly*ty)/B = A=pB 
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DISCUSSION 


In the case of finite geometries, the absolute neutron density at any 
point in space at any instance, is larger on the two group approximation 
than on the age approximation; for, on the age theory, the non-leakage 
term for fast neutrons is exp. [— tT w,”] whereas on the two group theory 
it is (1+ 7 w,*)-". However, as the size of the moderator decreases, 


exp. (— 7 9,”) approaches (1+ t9w,”)-! and hence there is better agreement 
between the two theories. 


In the case of ordinary water, the age and two group approximations 
have been applied for infinite geometries. The agreement of these results 
has been compared with experimental data by Ramanna et al.” 


We are grateful to Dr. R. Ramanna for suggesting the problem and 
for helpful discussions throughout the course of the work. We also thank 
Mr. Udgonkar for many helpful discussions. 


REFERENCES 


1. Gallone and Salvetti .. “Symbolical calculation methods applicable to neutron multi- 
plication,” JJ Nuovo Cimento, 7, 482. 
“Time-dependent neutron sources in multiplying media, II,” 
ibid., 8, 960. 
2. Ramanna et al. .. On the determination of diffusion and slowing down constants 
of ordinary water and beryllium oxide using a pulsed 


neutron source, Paper read at the Geneva Conference 
(to be published). 














Thermal Neutron Distribution in Non-Multiplying Media 187 
APPENDIX A 


SYMBOLIC SOLUTION OF THE TIME-DEPENDENT ‘ AGE” EQUATION 
Consider the balance of neutrons of energy between E and E + dE. 
We get 


Dvn (r, t, E)dE + 2, X(t, DJ dE = ~ (r,t, B) dE, (1) 
3E ry, 


where 


n(r, t, E) dE is the density of neutrons of energy between E and E + dE 
at the position r and at time ¢. 


X (r, t, E) dE is the slowing down density of neutrons at energy E at the 
position r and at time ¢ 


and 


D is the diffusion coefficient of the medium for neutrons of energy E. 
The following relations hold 


A 
n(r, t, E) dE= 2B X (r, t, E) dE; D = »; v/3. (2) 


where 
A;= scattering mean free path for neutrons of energy E in the medium. 
A¢= transport mean free path. 
v = velocity of neutrons of energy E. 
€ = average logarithmic loss of energy. 
Substituting in (1) 
AsAt 2 OX O71 As ox 3 
38 VY * TB BE os wie 


Introduce a new variable 7 called age defined by 


Eo 
To = 32 f As (E) Xe (B) (4) 


from (3) 
AsAt AsAt dX i As &X 


eS eee a 2. bow 
3 VY" — 328 3y, ~ SE 2B 
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or 


Vx (r, t, T) ~ me [x (r, t, 7) = Dm) - [x (r, t, T))]- (5) 


This is the time-dependent age equation. 
Let x (r, t, 0) - do (1, t) 


where ¢,)(r, ¢) can be expanded into orthonormal eigenfunctions under 
appropriate boundary conditions as previously described in the text. 


Then the solution is 


X (rr, t, 7) = exp i fen] V* -| f or i 3 o(r, 2) 


from (4) we can get 


t 
%= J ’D dt where ty is the slowing down time for a neutron of age 75. 
0 


To 1 ty T, 
Then f D(m) dt, = f dt = t;; f dry = 7%. 


So, X (x, t, 7) = exp [ TV" — ty si Go (F, ¢) 
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APPENDIX B 


EVALUATION OF THE INTEGRAL I (y, v) 


( e* sin kr 
Li )= f ka dk 
=f kev sin kr | Fema "kd dk 


e* dx f kere 1K? sin kr dk 


0 


I 
Oon8 


= AF dxf Ke mv"O J, (ker) dk 


F 7 [(u-v?2)) 


_ va or dines 
. fe Gx + we & 


x=0 t=r[2V/p 


Let t=r/2(u+ v*x)t; when 


Since x varies from 0 to oo, x is positive. Similarly ¢ is positive since it 
varies from 0 to r/2/u, r and +/p» being positive as they are physical 
constants. 


y- r[2V'p . 
I(u, v) = et” = f exp [- (aa +) dt 
: rl2/e rl2V a 
2 
~ for |er few [- (G- J ]arer f 
0 0 


exp [- (7 + +) ] dt | 
ae ye eulv® [e"'” I_+ et*T,] 


where 
ri2V 


u= f aol - (+4) ]4 
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Consider the integral I 


Let (5, —t ) = u. 


pa UE VP Ht 2r/v 
- 5 ; 





If u is negative, then we can define u’ = (t — r/2tv) 
so that 


t= 





uo + v/u® + 2r/v 
2 


since ¢ is always + ve and since /y? + 2r/v> u and Vu’? + 2r/v> wu’ 
in both cases, the positive sign must be used. 


The limits of integration are given by 


t=0—>u= oo and t= ~- 


rl2/ bh 


L= foo [ (S.-J) 


( - 57) (“ pe vi) 


7a) 
v 2 u 
ew udu —: girly ’ Pamela 


Vue + 2rJv ° Vay ole 


co 


Vu? + 2r/v 
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a = [1— en (MH - xi) | ais _ e2tiy [1 — ef (MM + mn) 


Similarly it can be shown that 


w= ere [1 — an — 85 )]— Ete ME + oh) 


16) = Pew ler + eb] 


=7 2 eel ete { 1— erf (Ve — x) 


—e {i-eg (¥ + x7) f 





SYNTHETIC EXPERIMENTS IN THE BENZOPYRONE 
SERIES 


Part LV. A New Synthesis of 7-Hydroxy-Chromeno-(3’ : 4’: 2: 3)-Chromone 


By A. C. MEHTA AND T. R. SgesHADRI, F.A.Sc. 
(Department of Chemistry, University of Delhi, Delhi) 


Received October 1, 1955 


IN an earlier publication! a general method for the synthesis of the four 
ringed chromeno-chromone unit of the rotenoid skeleton was described 
based on considerations of possible biogenesis from the related 2-methyl 
isoflavone. There were a number of stages in this synthesis and the yields 
were poor. Successful attempts have now been made to simplify the 
synthesis of such compounds in order to make them available in good yields 
for further synthetic work. 


The new method is based on the earlier findings of La Forge,? Baker and 
co-workers,? and Mehta and Seshadri* that acid chlorides can be made to 
condense with desoxy benzoins in presence of pyridine to yield isoflavones. 
Ethoxy acetyl chloride® has now been used as the appropriate reagent. A 
pilot experiment using w-phenyl resacetophenone® gave a good yield of 


2-ethoxymethyl-7-hydroxy isoflavone (I) which on treatment with hydro- 
bromic acid-glacial acetic acid mixture is converted into the 2-bromo- 
methyl derivative which has been reported earlier.? When 2: 4-dihydroxy-2’- 
methoxy phenyl benzyl ketone’ is used, 2-bromo-methyl-7: 2'-dihydroxy 
isoflavone (IIb) is obtained through the intermediate (IIa). It has been 


_— Hg YY Ve 
Who AK 
(a) R=OC,Hs; R'=OCH, 


He 
6) R=Br; R’=HO 


( ( os 
NANG Nee 


Ill 
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found during this work that treatment of 2-ethoxy methyl isoflavones with 
hydrobromic acid at higher temperature results in considerable resin forma- 
tion and good results are obtained by employing mild conditions as de- 
scribed in the experimental part. The bromo-methyl-dihydroxy isoflavone 
(11 b) readily eliminates a molecule of hydrogen bromide in the presence of 
anhydrous potassium carbonate in dry acetone giving rise to 7-hydroxy- 
chromeno-(3’: 4’: 2: 3)-chromone (III). 


EXPERIMENTAL 
2-Ethoxymethyl-7-hydroxy isoflavone (I) 


w-Phenyl resacetophenone (5 g.) in dry pyridine (30c.c.) was treated 
with ethoxy acetyl chloride (10c.c.) at 0°. The product of condensation 
was worked up after 24 hours in the manner reported earlier* for acetyl 
chloride. The intermediate diketone could not be obtained crystalline and 
pure. Hence it was refluxed with 10% aqueous sodium carbonate (200 c.c.) 
for 2 hours. The solution was cooled, filtered, and the filtrate acidified 
with dilute hydrochloric acid. The sticky solid that separated out was 
extracted with ether and the extract dried over anhydrous magnesium 
sulphate. The residue left after the removal of the solvent crystallised from 
a mixture of ether and light petroleum and finally from ether as colourless 
rectangular prisms, m.p. 146-47° (3 g.) (Found: C, 72-9; H, 5-2; C,;H,.O, 
requires C, 73-0 ; H, 5-4%). The mother liquor on concentration yielded 
the original ketone (1-2 g.). 


2-Bromo-methyl-7-hydroxy isoflavone 


The above isoflavone (I) (1 g.) was dissolved in glacial acetic acid (10 c.c.) 
and hydrobromic acid in glacial acetic acid (50%; 10 c.c.) added. After 
keeping at room temperature for 12 hours a further quantity of hydrobromic 
acid (50%; 10c.c.) was added and the mixture refluxed on a water-bath for 
45 minutes. Ice-cold water was added and the solid that separated out was 
filtered and dried. It crystallised from a mixture of ethyl acetate and light 
petroleum as tiny rectangular prisms, m.p. 246-47° (0-8 g.). Mixed m.p. 
with an authentic sample’ was undepressed. 


2-Ethoxymethyl-7-hydroxy-2'-methoxy isoflavone (II a) 


2: 4-Dihydroxy-2’-methoxy phenyl benzyl ketone (5 g.) in dry pyridine 
(30 c.c.) was condensed with ethoxy acetyl chloride (10c.c.) at 0° and the 
product worked up as in the previous case. 2-Ethoxymethyl isoflavone 
(II a) crystallised from a mixture of ether and light petroleum and finally 
from ether as colourless rhombohedral prisms, m.p. 168-70° (3-2 g.) (Found: 
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C, 69-3; H, 5-7; C,)H,sO; requires C, 69-9; H, 5-5%). The mother liquor 
on concentration gave the original ketone (1 g.). 


2-Bromo-methyl-7 : 2'-dihydroxy isoflavone (II b) 


The above isoflavone (II a) (1 g.) was treated with glacial acetic acid 
and hydrobromic acid as done in the simpler case. The product crystallised 
rom a mixture of ethyl acetate and light petroleum and finally from ethy] 
acetate as colourless stout rectangular prisms, m.p. 213-14° (decomp.) 
(0-8 g.) (Found: C, 55-4; H, 3-5; C,gsH,,O,Br requires C, 55-3; H, 3-2%). 


7-Hydroxy-chromeno-(3’ : 4’: 2: 3)-chromone (IIT) 


The above isoflavone (IIb) (0-3 g.) was dissolved in dry acetone 
(100 c.c.), freshly ignited potassium carbonate (3 g.) added and the mixture 
refluxed for 4 hours. More potassium carbonate (3 g.) was now added and 
the refluxing continued for 6 hours more. Acetone was distilled off, the 
residue dissolved in water and acidified with dilute hydrochloric acid. The 
solid product crystallised from alcohol as colourless needles, m.p. 240-42° 
(decomp.) (0-2 g.). Mixed m.p. with an authentic sample! was undepressed. 
The acetate obtained by the acetic anhydride-pyridine method also agreed 
in its melting with the authentic sample. 


SUMMARY 


By the use of ethoxy acetyl chloride and pyridine at 0°, 7-hydroxy- 
chromeno-(3’ : 4’: 2: 3)-chromone has been synthesised much more conveni- 
ently than by the earlier methods reported in the literature. 
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1. INTRODUCTION 


UsING a method based on series expansions in Mathieu functions, we have 
shown the following property concerning the intensities of the lines of the 
spectra caused by the diffraction of light by supersonic waves at normal 
incidence of the light': if only the zero and first orders are visible, the values 
of the intensities of the orders + 1 and — 1 for standing supersonic waves 
are half the values of the intensities of the same orders in the case of pro- 
gressive sound waves. ~ 


A confirmation of this property was found in the experiments of 
Sanders,2, which give a systematic comparison between the diffraction 
spectra caused by progressive and standing supersonic waves. 


The purpose of the present paper is to deduce this property from the 
generalised theory of Raman and Nath.*® 


2. THE RESULTS OF NAGENDRA NATH FOR PROGRESSIVE SOUND WAVES* 


According to the generalised theory of Raman and Nath the amplitudes 
of the diffracted light waves satisfy the following system of difference 
differential equations 


2 a — (¢na— dnis) = in*pdy (n = — 00,...., + 00), (1) 


with the boundary conditions 
$y (0) == 1, dn (0) =90 (n £90); (2) 
where 
__ Qampz 


A 5] 
z= the width of the sound field along the direction of the incident light, 


2 
, HpgA*? : 
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A = the wave-length of the incident light in vacuum, 
A*= the wave-length of the supersonic waves in the considered medium, 
/o = the refractive index of the undisturbed medium, 


# =the maximum variation of the refractive index. 


Nagendra Nath has further shown,‘ that if the parameter p> 1, only 
a small number of orders will appear. In the case that only zero and first 
orders are visible, the system (1) becomes, 


ee th=0 | Gu=— Fy) (3.4) 


2 oe — $o = iphy (3 5) 


and its solution satisfying the boundary conditions (2) is 


$o = [1— we 3 sin? fe + 8) g x 


P 24 8) 
col! (fae (Geto CEM] oe 
2 .. ((p? + 8) , ) _€ 
b= Grpapsin [Po ey exw | gy “- 
The first order intensity is then 


Oaee.. ; 2 {(p? + 8) 4 - oP ¢ 
h= +73 ) 4 | ~ oa sin 4° (5) 


3. THe RESULTS FOR STANDING SOUND WAVES 


In the case of standing sound waves the system determining the ampli- 
tudes of the different diffracted orders is® 


2 ar — sin € (¢?n1—$n+r) = in*pdn (n = —oo,...., + coc) (6) 


with 
e« = 2rv*t, (6 a) 


v* = the frequency of the sound waves. 


The boundary conditions still are 


$y (0) = 1, $0) = 9. (7) 
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If orders higher than the first are to be neglected, the system (6) becomes 
explicitly, 


d 
apt ds sin «e = 0 (8 a) 
5 4 ; ; 
2 — do Sin € = ipdy. (8 db) 
Making the substitution 
C= sin e (9 a) 
= p/sin €, (9 b) 
we can write this system as 
d 
- + $1 =0 (10 a) 
- ; 
a — $9 = iogy, (10 5) 


which has the same form as the system (3 a, b) and, of course, the same 
solution, if one respectively replaces ¢ and p by ¢ and o. Formula (4 5) 
giving the first order amplitude thus becomes, 


e (o? + 8) iog 
¢, = arp sin { z= ches 


_  2sme (p? + 8 sin? et ipt 
~ (p2 + 8 sin? e)t sin 4 é} wtibs 


or taking into account, that p > 1 
2sin iot e* — ete . ipt 
oy © 4s in (5 é)e% - —— sin (fe ) rs (11) 


This wave is thus splitted in two waves, having frequency-changes of the 
amount +v* and —v*. The amplitude of the wave undergoing a frequency- 
change +v* is Ay = (1/ip) sin (p€/4) exp (ipé/4), while the amplitude of the 
wave having a frequency-change —v* is A= (— 1/ip) sin (p§/4) exp (ipé/4). 
So the total intensity of the first order becomes 


iy () = RRs + AAM= 4 sin? (12) 


The property stated in the introduction now follows immediately from 
comparison of the formule (5) and (12). 
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4. REMARKS 


(1) If sin « = 0(t = m/2v*), it is seen from the system (8a, 5), that 
¢) = constant = 1 and ¢,= constant x exp (ipé/2) = 0, this last constant 
being zero according to the boundary conditions (7). 


(2) The formula (12) may also be obtained in another way, which does 
not take account of the spectral character of the different orders. From 
the formula (5) the first order intensity in the case of standing sound waves 
may be written 
(o? + 8) 4sinte 2 {(p? + 8 sin? «) 

$y (tee 


= Sl 


oS . 
h(@) = o? +8 an? p* + 8 sin?e 


4sin’e . 
wv SS sinth € = (p>). 
This intensity is depending periodically on time, so that only the mean 
intensity will be observed experimentally. This mean value is given by 


hig@=- hi = ‘ sin? (6), f sin? aver dt =~, sin? ¢ 


(7*= the period of the sound waves), which is the value given by (12). 
5. SUMMARY 


It has been shown that the following property of the spectrum of the 
diffraction of light by supersonic waves, in which orders higher than the 
first are missing, can easily be deduced from the generalised theory of Raman 
and Nath if p> 1: the intensity of the first orders in the case of standing 
sound waves has half the value of the intensity of the same orders in the case 
of progressive sound waves. 


The author wishes to express his thanks to Professor Dr. M. Nuyens 
for the valuable discussions he had with him. 
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TAMARIND seed jellose is a polysaccharide* which is composed of galactose, 
xylose and glucose in the ratio of 1: 2:3. 2 The repeating unit of the macro- 
molecule has been recently shown to have the following structure, containing 
16 free hydroxyl groups.* 
X 
1 G = Glucose 
2 Ga = Galactose 
Ga xX X = Xylose 
I 1 
6 6 


The acetylation of the substance has been carried out under the same condi- 
tions as those adopted for starch,* and the results are reported in this paper. 


Jellose, like starch, undergoes progressive acetylation, depending upon the 
acetylating agent and the conditions of the experiment. If one hydroxyl 
in every one of the six saccharide rings of the hexasaccharide unit undergoes 
acetylation, a hexaacetate results corresponding to the monoacetyl derivative 
of starch. If two of the hydroxyls in every one of the six saccharide rings are 
acetylated, the product is a duodecaacetate corresponding to the diacetate of 
starch. If all the hydroxyls are acetylated, a hexadecaacetate forms, which 
corresponds to the triacetate of starch. Corresponding to these stages, the 
hexaacetyl, duodecaacetyl and hexadecaacetyl derivatives of jellose have 


* In a recent note (Curr. Sci., 1955, 24, 235) Savur has claimed to have isolated from tamarind 
kernel powder three distinct types of polysaccharide fractions, differing in their solubility and 
power of gelation. The author appears to suggest that the above fractions are different chemical 
entities. However, it is very likely that they may be the same polysaccharide existing in different 
degrees of polymerisation. The material used in the present investigation as also that employed 
for the earlier constitutional studies (Rao, P. S. and White, E. V., J. Am. Chem. Soc., 1953, 75, 2617) 
was prepared and purified according to the procedure described earlier (Rao, P. S., Ghose, T. P. 
and Krishna, S., J. Sci. and Ind. Res., 1946, 4, 705). It did not show any heterogeneity on repeated 
precipitations from dilute aqueous solutions by means of alcohol, the purified product having a 
specific rotation of 71-4° in 0-25 per-cent. aqueous solution and a pentosan content of 30-5 
per cent. (Rao, P. S. and Krishna, S., Curr. Sci., 1946, 15, 133). 
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been obtained. It may, however, be noted that the reaction need not in- 
variably take place in the stepwise manner indicated; in fact, as in the case 
of starch, intermediate acetyl derivatives are also formed. 


It may be mentioned that the polysaccharide undergoes acetylation 
even with boiling glacial acetic acid in the presence of anhydrous sodium 


acetate (cf. starch®) and the extent of acetylation depends on the duration 
of the treatment. 


The acetyl derivatives are as a class insoluble in water and, some, even 
in organic solvents. The solubility appears to depend not only on the extent 
of acetylation but also on the degree of depolymerisation that the molecule 
undergoes during the chemical treatment. The insoluble products can be 
rendered soluble in organic solvents by treatment with a hot mixture of 
sulphuric and glacial acetic acids. In appropriate solvents some of the acetyl 
derivatives yield fairly strong, flexible, glossy and transparent films which 
adhere to glass, wood and metal surfaces. From this point of view the 
product obtained by the action of boiling acetic acid is particularly good. 


Further, some of the acetates, on account of their wide melting-range, may 
be useful as thermoplastic resins. 


EXPERIMENTA L 


Action of glacial acetic acid—Dry jellose (10 g.) was refluxed with 200 c.c, 
of glacial acetic acid and one gram of anhydrous sodium acetate for varying 
periods. After each treatment the mixture was poured on crushed ice and 
the resulting solid was filtered, boiled with small quantities of fresh water to 
remove any unreacted jellose, cooled, filtered, washed with alcohol and 


finally dried in vacuum. The properties of the different products are given 
in the following table :— 








TABLE I 
Length of P 
soontanant Acetyl value Properties 
12 hours | 12-39% | Insoluble in water and all the ordinary organic 
| solvents. 


24 hours 16-58% Behaves similar to the previous case but swells in 


glacial acetic acid. 


36 hours 19-98% | Soluble in alcohol, pyridine and dilute acetic acid but 


not in water, acetone and chloroform. 
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In the first two cases (boiling for 12 and 24 hours) the yields were not high, 
because the products were contaminated with considerable quantities of 
unreacted jellose. In the third case, however, the yield was 8 g. from 10 g. of 
the jellose. In this treatment the reaction product went into solution in the 
acetic acid and did not separate out on simple dilution, but only on satura- 
tion with sodium chloride. Its solutions on slow evaporation left fairly strong, 
pliable and transparent films. It softened at 270° and became dark at higher 
temperatures, but did not melt below 320°. 


Jellose hexaacetate.—The jellose was initially gelatinized by treatment 
with formic acid. Dry jellose (10 g.) was gradually added to 100c.c. of 90 
per cent. formic acid with continuous stirring. The substance immediately 
swelled, forming a thick viscous solution. An analysis of the isolated pro- 
duct showed that the formic acid reacted with the jellose, converting the 
latter into the hexaformyl derivative [(a))*° in formic acid solution 
(c = 0-3%) = + 46°3°; softens at 260°, swells enormously on further 
heating but does not melt below 320°]. Continuing the stirring, acetic 
anhydride (100c.c.) was added dropwise. Since carbon monoxide was 
evolved as a result of the decomposition of formic acid, the operation was 
carried out in a fume chamber. During the reaction the temperature rose 
gradually and had to be controlled by means of crushed ice. When the 
addition of acetic anhydride was over, the reaction mixture was cooled and 
treated with ice-cold water. On vigorous stirring a faintly brownish solid 
separated. The yield was 12g. It was insoluble in alcohol, acetone, 
chloroform, pyridine, etc., and was moderately soluble in formic acid, in 
which [a]p%° was + 59-4° (c= 1%). It softened at 230° and melted with 
decomposition at 280°. This corresponded to the hexaacetyl derivative of the 
jellose (Found :—COCHs, 22-4; [C3qH59O2s (COCHs)g]n requires—COCH,, 
22-2%). 


In the above experiment smaller quantities of formic acid could be used 
for swelling the jellose but the difficulty was as regards the stirring. Even 
in 10 per cent. concentration, the jellose formed a thick dough with the acid. 
However, when the formic acid was sprayed on the jellose triturating in a 
mortar, the swelling was quite rapid and efficient. By this means the weight of 
formic acid necessary for swelling could be reduced to an amount equal to 
the weight of jellose. 


Jellose duodecaacetate-—The experiment was carried out as in the pre- 
vious case but with the difference that acetic anhydride containing 1 per cent. 
of concentrated sulphuric acid was used instead of acetic anhydride alone, 
and the product was duodecaacetate. The isolation and purification were done 
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as in the previous case (Found:—COCHs, 36-0; [C3,4HssO2g (COCHS),9}n 
requires—COCHs,, 36°4%). The yield was 12-5 g. from 10 g. of the jellose. 
It was soluble in acetone, acetic acid and pyridine and was insoluble in chloro- 
form and alcohol. In pyridine solution (0-75%), [a]p%° was + 51-6°. 
Softening at 240°, it did not melt below 320°. 


Jellose hexadecaacetate——If in the above experiment the reaction 
mixture was heated at 95°C. for two hours with continuous stirring, the 
completely acetylated product was obtained. During the first hour the 
temperature was slowly raised to 95° and was subsequently maintained 
thereat for the rest of the period. On dilution with ice-cold water and stirring, 
the acetyl derivative separated out. The yield was 16-6 g. from 10 g. of the 
jellose. The product was soluble in acetic acid, acetone, chloroform and 
pyridine but not in alcohol, and exhibited a specific rotation of + 28-5° in 
pyridine solution (c = 0-5%). It softened at about 240°, assumed a dough- 
like structure at 280° and did not melt below 320°. It was largely the completely 


acetylated derivative (Found:—COCH;, 44:3; [C34H4Oes (COCHS),6]n 
requires—COCHg, 43-49%). 


Conversion of insoluble into soluble acetates—The insoluble acetate 
(10 g.) was suspended in 100 c.c. of glacial acetic acid containing 1-5 per cent. 
of concentrated sulphuric acid, and the suspension was heated on an oil- 
bath. The temperature was slowly raised to 95° during the course of an 
hour and the heating was continued for another hour at the same temperature. 
The mixture was then diluted with ice-cold water, whereupon the “ soluble ” 
acetate separated (yield: 9-5 g.). 


Experiments on the utilisation of the acetates——Since the melting range 
of the acetates is fairly wide, covering nearly 60° from the softening point to 
the formation of the melt, the substances, it was felt, might be useful as 
thermoplastic materials. However, preliminary experiments on moulding 
of the esters at temperatures just a little above the softening points indicated 
that the casts were rather brittle, though the surface was glossy. Addition 
of plasticisers like camphor, phthalic anhydride and esters of phthalic acid in 
varying amounts improved the properties of the casts to different degrees, 
Phthalic anhydride and esters of phthalic acid (25 per cent. addition) gave 
the best results but even then the casts were not sufficiently strong and 
tough. In these experiments the duodeca ester gave much better results 

than the other acetates. The addition of cellulose acetate by about 25 
' per cent. improved the strength of the jellose acetate moulds (cf. starch 
esters‘). 
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The films formed by the evaporation of the solutions of the hexadeca- 
acetate, though transparent and glossy, were very brittle. Addition of 
plasticisers did not improve the properties appreciably. On the other hand 
the films formed from the duodeca ester were strong and tough and plasticisers 
(the anhydride and the esters of phthalic acid), when added in about 25 
per cent. further improved the properties. Being transparent, glossy and 
water-resistant, the films may find application in the preparation of water- 
proof covers. In this connection special mention may be made of the lower 
acetylated products obtained by the treatment of the jellose with boiling 
glacial acetic acid for 36 hours. Their films from pyridine solution were 
pliable, fairly strong and transparent. 


When applied to wood in pyridine solution the coatings of plasticised 
acetates resembled varnishes. These coatings adhere to glass and metals 
as well, and on drying, particularly at higher temperatures, they become 
flexible and impermeable to moisture. A suitable concentration in these 
experiments was found to be 5 per cent. 


SUMMARY 


The hexa, duodeca and hexadeca acetyl derivatives of tamarind seed 
jellose have been prepared by treatment with acetic anhydride under different 
conditions. Boiling glacial acetic acid also effects partial acetylation. 


Some of the acetyl derivatives in appropriate solvents give fairly strong, 
flexible, glossy and transparent films, which adhere to glass, metallic and 
also wooden surfaces. They may also be useful as thermoplastic resins on 
account of their wide melting range. 
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ABSTRACT 


The lunar semidiurnal variation in the meson intensity at Kodai- 
kanal has been determined by the Chapman-Miller method of analysis 
from bihourly records of the cosmic ray meson intensity registered by 
three-fold coincidences for two consecutive groups of data, each extending 
over a period of 18 months. The former period gives a result about two 
to three times the probable error and can be taken to be significant. 


While the correlations between the solar diurnal variation of meson 
intensity and barometric pressure are high and negative, the correspond- 
ing correlations between the lunar semidiurnal variations of meson 
intensity and barometric pressure are positive. The results obtained at 
Kodaikanal agree with those of Duperier obtained at London. 


But the latter period does not give any significant result outside the 
probable error. Hence further investigation has to be continued. 


I. INTRODUCTION 


THE lunar tide produced in the atmosphere, though small compared to its 
counterpart solar tide, is nevertheless important. It is revealed in the 
lunar semi-diurnal variations of the barometric pressure,! of air temperature,? 
of the height of the E-layer,* of wind components,* of geomagnetic elements*® 
and of the cosmic ray meson intensity registered at ground level.6 In addi- 
tion to producing a mass absorption effect in cosmic rays, it is to be expected 
that the pressure oscillation would result in a vertical shift of the isobaric 
levels. If this shift is appreciable at the height where the bulk of the cosmic 
ray mesons are supposed to originate, there would result an increase or 
decrease of the meson intensity at the ground level, depending on whether 
the level of formation goes down or up. 


A study of the influence of the lunar tide in the atmosphere on the 
meson intensity of cosmic rays, is important not only for understanding the 
atmospheric oscillations at various levels, but also for checking our current 
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theories on the transition effects of cosmic rays in our atmosphere. So far, 
the only determination of the lunar semidiurnal variation L, in meson 
intensity is the one reported by Duperier* from data collected at London. 
He has reported an amplitude twice the probable error, and while this is 
indicative of the existence of L,, the evidence is not conclusive. 


Duperier is reported to have given further results of the effect of 
atmospheric lunar tide on the cosmic ray intensity, at the Cracow conference 
in 1947. His later result is reported to have an amplitude nearly thrice 
the probable error. The positive correlation between the L, for meson 
intensity and L, for barometric pressure is also surprising, in contrast to the 
negative effect observed in the solar daily variation. Since at low latitudes, 
the lunar and solar daily variations of barometric pressure are much larger 
than at high latitudes, it is particularly appropriate to study their effects 
on meson intensity at an equatorial station like Kodaikanal. (Latitude 
10° 14’ N., Longitude 77° 28’ E., Geo-magnetic Latitude 1°N., Altitude 
7,688 ft. above sea-level). 

II. PROCEDURE 


The lunar daily variations are generally masked by irregular variations 
and by solar daily variations which are comparatively larger. However, 
the elegant method developed by Chapman and Miller’? enables the deter- 
mination of the effect of the lunar atmospheric tide on the meson intensity. 
The cosmic ray meson intensity has been collected at Kodaikanal during 
the period May 1952 to December 1953 by one of the authors (D. V.) and 
from January 1954 to June 1955 by the other (T. S. G. S.). Separate records 
have been obtained for vertical intensity M, of radiation penetrating 7 cm. 
of lead absorber and for intensity M, of radiation penetrating 15cm. of 
lead absorber. Both M, and Mg correspond principally to intensity of 
mesons, though the latter contains relatively more energetic particles than 
the former. The mean bihourly rate for meson intensity is about 4,000 
triple coincidences per bihour per telescope. 


Though hourly photographs were available only bihourly data have 
been utilised, because, this reduces the labour to half with hardly any loss 
in the accuracy of the results. The days are sorted out according to their 
» values, where » refers to the lunar phase. In the course of a synodic 
month, » decreases from 24 to0. » = 24 or 0 is new moon, pu = 18 is first 
quarter, » = 12 is full moon and u = 6 is the last quarter. For the Chapman- 
Miller-method of analysis, only 12 groups are needed. Days of yp value r 
are combined with days of » value (r+ 12). Incomplete daily sequences 
have not been taken into consideration and only days for which there is no 
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break have been used. The daily sequence consists of thirteen bihourly 
terms, the thirteenth being the value for the initial bihour of the succeeding 
day. Each daily sequence is a unit and can be combined with other daily 
sequences. Summation of the sequences corresponding to the twelve 
groups of » is carried out. This method of using a thirteen term sequences 
for each day allows for a correction to be applied for the average non-cyclic 


change. The analysis has been done on the same lines as described by 
Tschu.*® 


A preliminary calculation with the data of one of the authors (D. V.) 
showed that when all the days are used, the probable error turns out to be 
very large. It is therefore necessary to select the days. This procedure 
is necessary as irregular large fluctuations tend to mask the small periodic 
variation under search as has been pointed out by Chapman. For this 
purpose, oc, the expected standard deviation in a bihourly estimate, is calcu- 
lated. o = 1/n, where n is the mean bihourly counts. The actual bihourly 
deviation, 4, from mean of each daily sequence is then compared with o 
to provide a criterion for sorting the daily sequences in terms of the degree 
of disturbance. 


The sequences are sorted out into three groups 4 +20, 4 + 3a, and 
A > So respectively. The major portion of the data falls within the second 
group. The number of days falling within the first group are unevenly 
distributed among the various » groups and hence is not useful for study. 
When the third group is taken, the probable error turns out to be large. 
Hence it was found that days for which 4 + 30 would be a good criterion 
for selection. 


In the process of calculating the lunar semidiurnal variation L,, the 
solar semidiurnal variation S, is also derived. This can be compared with the 
value derived directly from the 12 solar bihourly deviations from mean for 
sequences for all » values grouped together. 


For atmospheric pressure only days in which the difference range did 
not exceed 3 millibars have been taken into consideration. Although the 
lunar semidiurnal variation of the barometric pressure at Kodaikanal has 
already been calculated,® it has been redetermined for the period of experi- 
mentation, as a check on our method of calculation. 


An alternate method of determining the lunar semidiurnal variation 
L. in meson intensity and in the barometric pressure at Kodaikanal has also 
been attempted. The daily sequences are separated into 24, groups 
instead of 12-0 hour of a day for which p is 24 is nearly equivalent to 2 hours 
of a day for which p is 22 and is nearly equivalent to 4 hours of a day for 
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which » is 20 and so on. So the data corresponding to 0 hour of days 
of » = 24 are combined with those of 2 hours of days of » = 22 and then 
with those of 4 hours of days of » = 20 and so on. Similarly for all other 
bihours, addition is done to give the net summation of the group sequences 
according to lunar hours 0, 2, 4, etc. 


The the odd group, namely days of » values 23, 21, ...., 1, are taken 
and addition is carried out in the same way as pointed out in the previous 
para. But 2 hours of a day for which yp is 23 will be be equivalent to | hour 
of a day for which p is 24, 4 hours of a day of » = 23 will be equivalent 
to 3 hours of a day of » = 24 and so on. Hence from this group we get 
the net summation of the group sequences according to odd lunar hours 
1, 3, 5, etc. These 24 values are subjected to Fourier Analysis to determine 
the second harmonic. 


Ill. RESULTS 
A. Data from May 1952 to December 1953 


Table I gives the distribution of the days according to » groups, for 
meson intensity M, for 4 +20, 4 +30, and 4 +5oe; for M, for 4 + 2a, 
4 +30, and 4 +50; for (M,+M,) for 4 + 3c. 


TABLE I 


Distribution of days according to w values 








nN oF 7 r=) © ~ re) a S = a of = 
‘ be — — _ = fal =~ = N Nn N N on 
p values 2 ee & 8 3 & 2B B & oH a) Ee) oe 
oO -_ N of ~~ 9 © ~ @ oa So — BU 
= _ 
A>2e.. 6 10 12 16 14 14 12 15 ll 18 9 15 152 


Mi A>+3c.. 30 23 29 29 31 31 30 33 32 38 25 33-364 
A>5e.. 34 41 37 39 38 37 38 45 38 43 34 39° 463 





Ay>2e-- 12 14 13 12 12 17 19 21 12 15 19 15 18] 
Mz A>+>3c.. 26 27 28 29 29 33 33 37 32 35 32 26 367 
Arde 28 39 32 40 32 39 36 41 38 43 33 24 425 





(M;+Mz2) 56 50 57 58 60 64 63 70 64 73 57 59 = 731 
A>+>3e 





Table II gives the determination of 1,, the amplitude and A,, the phase 
of the semidiurnal variation L, for meson intensity M,, M, and (M, + M,) 
for the various cases and for pressure P. The probable errors are also given 
in the table. It can be seen that only for meson intensity M, for 4 + 3c, 
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TABLE II 
L; 
l, PE. As 

A>+20 0-07 0-07 —34° 34’ 

M, 4>3e 0-15 0-05 — 9 2’ 
4>5e 0-12 0-05 a 

4>20e 0-08 0-07 —54° 29’ 

M, 4>3c 0-04 0-04 —26° 35’ 
4>+5e 0-05 0-05 Ss 

(M, + M,) 4>+>30e 0-09 0-04 —14° 39’ 
P 40 12-8 31° 16’ 

microbars 





M, = Meson passing through 7cm. of lead absorber. 

M, = Meson passing through 15 cm. of lead absorber. 

P = Barometric pressure. 

4 = Degree of disturbance of the days chosen. 

P.E. = Probable error. 

I, = Amplitude of lunar second harmonic of meson, and pressure. 
Ag = Phase of lunar second harmonic of meson and pressure. 


M, for 4 + 5o and for (M, + M.) for 4 +30, and for pressure, are the 
results fairly significant. ; 


Fig. 1 shows on a lunar 12-hourly harmonic dial, the determination of 
L, for meson intensity M,, Ms, (M, + M,) and pressure P. The probable 
error circles are also indicated at the end of the vectors. 


Table III shows the value of solar semidiurnal variation S, derived 
during the process of calculating L,, along with the value derived directly 
by the usual method without separation of the days according to m» groups. 
S, refers to the solar semidiurnal amplitude and , to the phase. 
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Fic. 1. Lunar 12-hourly harmonic dial showing the lunar semidiurnal component 


meson intensity M,, M, and (M, + M,), and pressure P. 








TABLE III 
Ss 
Se be 
© ATG = 1-39% S20 —5% 
B 43% —68° 
S, = Amplitude of solar semidiurnal component of meson. 
be = Phase of solar semidiurne| component of meson. 
A = Derived during the course of the derivation of the lunar tide. 
B = Derived in the usual way, without separating the data according to lunar phase angle. 


for 


The 24 values obtained for M, and (M,+M,) for 4 +3c and for 
pressure P are smoothened by taking moving average of three consecutive 
values. These smoothened daily variation curves according to lunar hours 
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are shown in Fig. 2. The positive correlation between the lunar daily 
variations of meson intensity and pressure is clearly visible. 
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Fic. 2. Lunar daily variation curves of meson intensity M,,(M, + M,), and pressure P 
after taking the moving average of three consecutive hourly values. 

As a check on the Chapman-Miller method L, has been determined 

for meson intensity M, for 4 + 30 and pressure P by the alternate method. 

This is shown in Table IV. There is fair agreement between the two results. 


Table V compares the study at Kodaikanal during the period May 1952 
to December 1953, with Duperier’s study at London during the period 
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TABLE IV 
L, 
(Alternate Method) 
I, Az 
M, 0-14% —21° 
P 47 39° 
microbars 
TABLE V 
Greenwich Kodaikanal 
Meson Sy» 0-14% 0-43% 
fo —6° —68° 
Press. Sz 0:22% 0-74% 
po 144° 32> 
Meson /, 0-023% 0-15% for M,, 0-09% for (M, + M,) 
re 160° —9° for M,, —15° for (M, +> M,) 
Press. /, 0-009% 0-03% 
A, 114° 39° 
Meson S,/l,_ .. 6 3 for M,, 5 for (M, + M,) 
Pus. Sls = «- 24 25 
Yury 0-69 0-67 for M,, 0-59 for (M, + M,) 
Brae) 1-76% per mm. 3-35% per mm. for M,, 
1-77% per mm. for (M, + M,) 
YS MP) —0-87 —0-95 
Bsr, —0-55% per mm. —0-54% per mm. for (M, + M,) 
Sg = Amplitude of solar second harmonic. 
be = Phase of solar second harmonic. 
ly = Amplitude of lunar second harmonic. 
Ag = Phase of lunar second harmonic. 
‘Yuwup) = Lunar correlation coefficient between meson and pressure. 
Bue) = Lunar pressure coefficient. 
Ystur) = Solar correlation coefficient between meson and pressure. 
Bsiue») | = Solar pressure coefficient. 
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TABLE VI 
L, 
(Chapman-Miller Method) 
Meson l, P.E. rz 
% Yo 
M, 4>3c 0-04 0-04 a — 11° 
M, A>+>3e 0-05 0-04 am— 5° 
(M, + M,) 4>3e 0:04 0-04 a— T° 
TABLE VII 
L, 
(Alternate Method) 
Meson lps rz 
M, 4>+>3c 0-04 a — 25° 
M, 4>3e 0-05 a— 5° 
(M, + M,) 4>30e 0:04 aw — 13° 





May 1941 to April 1944. It is satisfactory that there is close similarity 
between the observed facts of the two places. 


B. Data from January 1954 to June 1955 


Only days for which 4 did not exceed 3 are chosen. Table VI gives 
the L, determination for M,, Mz, and (M, + M,) by the Chapman-Miller 
method, and Table VII gives the same by the alternate method. There is 
fair agreement between the two. 





IV. DISCUSSION OF RESULTS AND CONCLUSION 
A. Data from May 1952 to December 1953 


For barometric pressure, the accurate determination of L, at Kodai- 
kanal by S. K. Pramanik® from data covering a period of six years gives the 
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amplitude 1, = 52 microbars and phase A, = 25°. The agreement of these 
with the values determined for the limited period of 18 months during which 


the cosmic ray experiment at Kodaikanal was conducted by one of the 
authors (D. V.) is quite good. 


The agreement between the values of the solar semidiurnal variation S, 
determined directly from bihourly deviations according to solar hours and 
those determined during the process of calculation of L, from data grouped 


according to lunar age as seen from Table III lends confidence to the L, 
determination. 


From the comparison of the results obtained at Kodaikanal and at 
London, the following points emerge :— 


(a) The 1, and s, for pressure increase by about 3-3 times from latitude 
52°N to latitude 10° N. In meson intensity the increases are of similar 
magnitude, being 3-9 for 1, and 3-1 for s,, for (M, + Mg) with 4 +3. 
For M, with 4 + 3c, the increase in 1, is however larger. 


(6) Both at London and Kodaikanal the correlation between the S, 
for meson intensity and S, for barometric pressure are high and negative. 
However the corresponding correlation between the L, for meson intensity 
and that for barometric pressure are positive, being + 0-69 at London and 
+ 0-67 at Kodaikanal. Thus we have the curious result of having a nega- 
tive barometric pressure coefficient for the solar tide S, and a positive baro- 
metric pressure coefficient for the lunar tide L,. In discussing his result, 
Duperier has neglected this important feature. He has moreover treated 
the amplitudes of the S, and L, variations without consideration of the 
striking change in phase in relation to pressure oscillation in the two cases. 


Therefore his remarks directly comparing the nature of S, and L, variations 
in meson are not valid. 


Duperier’s later result does not differ appreciably from his.earlier result, 
except that the latter is more significant. Therefore the conclusions drawn 
from his earlier result are valid for his later results also. 


B. Data from January 1954 to June 1955 


It is rather surprising that the analysis for this period has failed to give 
any significant results. All the amplitudes are of the same order as the 
probable error and as such no significant conclusions can be drawn from the 
same. It is difficult to understand why out of two consecutive periods, one 
should fail to give a significant result. It becomes therefore necessary to 
analyse further data before one can draw any definite conclusions about 
the effect of the lunar atmospheric tide. 
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